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We present a study of the Dirac eigenvalue spectrum near the region of the QCD phase transition. 
This study makes use of a sequence of ensembles with temperatures from 150 MeV to 200 MeV 
generated with 2+1 flavors of dynamical domain wall fermions (DWF) and the dislocation sup- 
pressing determinant ratio (DSDR) action on a 16^ x 8 lattice with an extent of 32 or 48 in the fifth 
dimension. All the simulations lie on a line of constant physics with 200 MeV pions. The DWF 
Dirac operator is normalized using the methods of Giusti and Luscher combined with those of 
Rome-Southampton collaboration, allowing a direct evaluation of the Banks-Casher relation. The 
relation between the resulting temperature-dependent Dirac eigenvalue spectrum and the possible 
restoration of f/(l)^ symmetry with increasing temperature is discussed. 
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1. Introduction 

The chiral phase transition is one of the most fundamental features of QCD. Lattice field 
theory has been applied successfully to the study of this interesting phenomena and the associated 
symmetries. While traditional lattice techniques measure the chiral observables in a straightforward 
manner, examining the low-lying part of the eigenvalue spectrum of the Dirac operator can provide 
unique insights into various aspects of the symmetry breaking or restoration that accompany the 
phase transition. 

For instance, the chiral condensate, the chiral order parameter for the transition, can be ex- 
pressed in terms of the eigenmode density via the following relation: 

/2m 
dAp(A) ^^^'l^^ ' ^ = ^^, (1-1) 

where p(A) is the spectral density of the Dirac operator and is the quark mass. When the chiral 
and infinite- volume limits are taken, one will obtain the well-known Banks-Casher relation [|l]], 

lim lim —{yiYi) = tt lim lim lim p(A). (1.2) 

In lattice calculation one may also examine the subtracted chiral condensate defined as, 

A/,, = (yiYi) (v>:vV^v) • (1-3) 

nig 

The subraction removes the ultraviolet divergent piece of the chiral condensate which is linear in 
quark mass. 

The Dirac eigenvalue spectrum can be utilized to study the anomalous U{1)a symmetry as 
well. A similar order parameter A;j_5, which is the difference between the pseudoscalar and scalar 
susceptibilities, can also be related to the eigenvalue spectrum. 

The expression above suggests that if there is a finite region above zero where the eigenvalue 
density vanishes (a gap), the U{1)a symmetry might be effectively restored. 

With chiral symmetry under good control, domain wall fermions (DWF) [§, ^ are an opti- 
mum tool for the exploration of the phase transition region. The residual chiral symmetry breaking 
(present in the DWF formalism for finite fifth-dimensional extent L^) is reflected in an additive 
correction to the bare quark mass (wires), which can be further suppressed by the adoption of the 
dislocation suppression determinant ratio (DSDR) action [m5|,|6|]. Despite some unphysical, mas- 
sive degrees of freedom from the extra fifth dimension, the low modes of the DWF Dirac operator 
should resemble an ordinary four-dimensional discretized version of QCD. Moreover, within the 
DWF formalism, the U{1)a symmetry is only broken by axial anomaly, rather than spoiled by 
lattice artifacts as with staggered fermions. 



2 



Dime Eigenvalue Spectrum at Finite Temperature Using DWF 



Zhongjie Lin 



2. Implementation Details 

We have collected eight ensembles near the phase transition region with 2 + 1 flavors of 
fermions. All the simulations have a 16^ x 8 space-time volume and a fifth dimension of Ls = 32 
or 48 and they all lie on a line of constant physics with rrij^ « 200MeV and kaon mass almost phys- 
ical [0]. Table I gives the basic parameters of these finite-temperature ensembles. The Ncfg column 
hsts the number of configurations for which the eigenvalues are calculated. Figure |I| shows the 
disconnected susceptibilities of various temperatures and indicates a critical temperature around 
160 MeV. 



T (MeV) 


/3 


Ls 








J3 




ji 


Wctg 




140 


1.633 


48 


0.00612 


-0.00136 


0.0519 


6.26(12) 


7.74(12) 


36(3) 






150 


1.671 


48 


0.00296 


0.00173 


0.0500 


6.32(29) 


6.10(29) 


41(2) 


340 


1.980(7) 


150 


1.671 


32 


0.00648 


-0.00189 


0.0464 


8.39(10) 


7.06(10) 


44(3) 


340 


1.905(6) 


160 


1.707 


32 


0.00377 


0.000551 


0.0449 


5.25(17) 


4.83(17) 


43(4) 


408 


1.725(8) 


170 


1.740 


32 


0.00209 


0.00175 


0.0427 


4.03(18) 


2.78(18) 


35(5) 


239 


1.631(11) 


180 


1.771 


32 


0.00132 


0.00232 


0.0403 


3.16(15) 


1.56(15) 


25(4) 


246 


1.476(4) 


190 


1.801 


32 


0.00076 


0.00258 


0.0379 


2.44(9) 


0.71(9) 


11(4) 


374 


1.439(3) 


200 


1.829 


32 


0.00046 


0.00265 


0.0357 


2.19(8) 


0.47(8) 


10(3) 


710 


1.365(3) 





1.750 


32 


0.00188 


0.00300 


0.0370 




252 


1.5685(5) 



Table 1: Summary of ensembles and the renormalization factors for the eigenvalue density. 
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Figure 1: Disconnected susceptibilites. Figure 2: Renormalization factors for the j8 = 

1.750, 16^ X 16 ensemble. 



For comparison, results from a zero-temperature ensemble with a volume of 16^ x 16 are 
presented as well. In order to keep = 200 MeV as the temperature and j8 decreases and mpes 
grows we must either increase L., above 32 or use a negative input quark mass. As the second 
and third lines of Table |T] show, this first use of a negative DWF input quark mass was successful, 
resulting in no exceptional configurations and giving a consistent result for X/,disc- 

We used the KaLkreuter-Simma [||] method to calculate the lowest A^eig = 100 eigenvalues of 
the Hermitian version of the Dirac operator Dh = /?^7^Ddwf> where is the reflection operator 
in the fifth dimension. Note the eigenvalues measured here include the mass term and are denoted 
by A, to be distinguished from those in the equations 1.1, 1.2 and 1.4. 
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3. Eigenvalue spectrum and its renormalization 

Renormalization must be applied to the eigenvalue spectrum of the DWF operator Dh before 
any sensible comparison can be made with either the input DWF quark mass m f or the eigenvalue 
densities from other four-dimensional fermion formalisms (e.g. Wilson fermions). The method 
we have adopted is a generalization to DWF of that proposed by Giusti and Luscher which 
introduces into the Lagrangian a twisted mass term, 

■^im{x) = Y.q'ix) (7^Dv + m + /Ai7^T')^^ W. (3.1) 

Then Green's functions such as the six-point correlator of the charged pseudoscalar density operator 
= q' {x)y^T^q'' (x) can be expressed in terms of the eigenvalue density, 

f^3(M) = -Lfc(^l)^2:3fe)/t4fe)A:5(^4)/t6fe)^6:i(^6)) (3.2) 



(3.3) 



t((75D)2 + ^2)3j/ 

dAp(A)(A2 + ju2)-3, (3.4) 

in the notation of Giusti and Luscher [Q]. The charged pseudoscalar density is well defined in the 
continuum limit in a variety of regularization schemes which are related by multiplicative renor- 
malization factors. Equation 3.4 implies the same rules also will work for the eigenvalue density 
e.g.: 

P!i' = ^Pl,' =^ P'i^')=Zml,n'p(^) (3-5) 



7 

... ... ^^n. .... , 

where Zy„^,„i is the renormaUzation factor for the mass term treated symmetrically with Pj^,. 

With such inspiration, we can invent a five-dimensional analogue of the twisted-mass term 
Pjj'^^ 'ix) = £fio' ^i{x,s)"^t'^i'{x,Ls — \—s) and relate it to the usual pseudoscalar density, 

1 

W{x)l^¥{x) ~ ^ £ ^'{x,s)y'^'{x,L, -\-s), (3.6) 

where i/a(jc) is the four-dimensional operator while ^(x,s) is the five-dimensional field. ' This 
renormalization factor Ztw^,,,/ connects the five-dimensional eigenvalue density to a more conven- 
tional density normalized in a fashion consistent with the usual bare quark mass mj: 

p'"f{A"'f)=Z^J_^,„^p{A^''y), A'"f=Z,„^A^''^ . (3.7) 



Because of the equivalence of the two operators at long distance expressed by equation [3.6| , 
the renormalization factor can be obtained from the ratio of the two Green's functions: 



The explicit sum over the fifth (s) dimension is suppressed later on for simplicity if no confusion is caused. 
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where T,'s are the PauU matrices in the flavor space. Results from a Coulomb gauge fixed wall 
source are presented in Table |l|. 

An alternative approach is to examine the off-shell, three-point Green's functions evaluated 
in Landau gauge. This is very similar to the Rome-Southampton non-perturbative renormaliza- 



tion (NPR) technique (RI/MOM) QlOp. The renormalization factor is extracted from the ratio of 



amputated vertices for the five- and four-dimentional operators. 

7(M0M) . . ^ Tr(l,,,,e'(^^^^-P-^')vA,fem(0)/?5y^^/'(0m^(^i)) 
^,„APUP2) Tr(l,,,,.'-(/'^^^-/'.-v.)v,,(^2)¥/(0)fw(0)v7/'(^i)> ■ ^ ^ ^ 

To fully utilize the whole lattice, we use a series of fixed-momentum volume sources to calculate 
the propagators, which is defined as 

T^ix;p)=e'P-n4y,4^hx3- (3.10) 

We perform our calculation using both non-exceptional kinematics, where p\ = p\ = {p\ — p2f' , 
and exceptional kinematics, where p\ = p2. The results for the zero-temperature ensembles are 
presented in Table |l| and Figure ^. 

Both Equation and should give consistent results independent of temporal separation 
t and of p\ and p2 respectively. Unfortunately, the NPR calculation is not feasible for the finite 
temperature ensembles due to large fluctuations. Therefore we only present the NPR results for 
the 16^ X 16 emsemble in Figure ^ Figure ^ shows a discrepancy between the two kinematics 
which is positively related to the physical momenta. This contradicts our expectation but can be 

plausibly explained by appreciable finite lattice spacing errors {apY at large momentum. Because 

in) 

the quantity Zj^i^^^. involves the smallest momenta, we use it to renormalize the spectrum. 

Further renormalzation from the bare mj scheme to the conventional, continuum MS scheme 



can then be easily performed since this final step has already been studied in detail [11]. Technical 
details and updated results will be available in our upcoming paper JT2|]. 

Figure ^ displays the effects of the renormalization, which can be naively regarded as a rescal- 
ing of the axes. The orange vertical line denotes the smallest of the hundredth eigenvalue and the 
spectrum below that is supposed to be complete. The other vertical lines indicate the bare masses 
of the light and strange quarks. The horizontal lines are the chiral condensates divided by 7i, which 



should agree with the eigenvalue density at A = as predicted by Banks-Cahser relation |LJ. There 
are two significant features associated with the renormalization. First, the light quark mass now 
matches the likely zero-mode peak. Second, the Banks-Casher relation agrees better although it 
is still inaccurate at 30% level. We attribute the discrepancy to finite-volume and finite-mass ef- 
fects. Thus, no definitive conclusion can be drawn before studies on a larger lattice and the chiral 
extrapolation are performed. 

Figure ^ shows the renormalized eigenvalue spectra at various temperatures near the phase 
transition region. Although not in perfect agreement with the Banks-Casher relation, at lower tem- 
peratures such as 150 and 160 MeV the chiral condensates and the eigenvalue densities are different 
from zero, signaling spontaneous chiral symmetry breaking. Above 170 MeV, these two quanti- 
ties both start to vanish as expected for temperatures above the transition. However, it remains 
uncertain whether the slope of the eigenvalue density vanishes before 190 and 200 MeV, where a 



5 



Dime Eigenvalue Spectrum at Finite Temperature Using DWF 



Zhongjie Lin 




Figure 3: Dirac eigenvalue spectrum for the j3 = 1 .750, 16^ x 16, zero temperature ensemble. The density in 
the left-hand panel has not been renormalized while that on the right has been changed into the normalization 
scheme of the usual input DWF mass nif. The left-most vertical line locates the total bare quark mass, 
nif + nires which matches well with the small peak seen in the renormalized, right-hand spectrum. 



possible gap does begin to emerge, indicating an effective restoration of ?7(1)a symmetry. A small 
peak at the lower end at these temperatures suggests that the major contribution to U{1)a symme- 
try breaking may come from zero-modes, which are expected to go away as the volume increases. 
Therefore, we await a calculation on a larger lattice to confirm this conclusion. 



4. Conclusions 

With the chirally symmetric DWF framework, we are able to explore the chiral and U{1)a 
symmetries near the phase transition region. The successfully renormalized eigenvalue spectrum 



of the Dirac operator as well as the correlator measurements [13] suggest an effective restoration of 
the U{1)a symmetry at temperatures higher than T^. We are looking forward to a similar simulation 
with larger volume to confirm our findings. 

I very much appreciate the help and advice from members of HotQCD and my colleagues at 
Columbia University. This work was supported in part by U.S. DOE grant DE-FG02-92ER40699. 
The simulations were carried out on the BG/P machine at LLNL, the DOE- and RIKEN-funded 
QCDOC machines and NYBlue machine at BNL. 
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